X-Ray Magnetic Circular Dichroism at the K edge of Mn 3 GaC 



Manabu Takahashi 
Faculty of Engineering, Gunma University, Kiryu, Gunma 376-8515, Japan 



m 
O 
o 

(N 



(N 



c3 



i 

-a 
c 

o 
o 



> 

(N 

in 
o 
m 
O 

-i— > 
c3 



i 

c 

o 
o 



X 



Jun-ichi Igarashi 

Synchrotron Radiation Research Center, Japan Atomic Energy Research Institute, Mikazuki, Sayo, Hyogo 679-5148, Japan 

(Dated: February 2, 2008) 

We theoretically investigate the origin of the x-ray magnetic circular dichroism (XMCD) spectra 
at the K edges of Mn and Ga in the ferromagnetic phase of Mn3GaC on the basis of an ab initio 
calculation. Taking account of the spin-orbit interaction in the LDA scheme, we obtain the XMCD 
spectra in excellent agreement with the recent experiment. We have analyzed the origin of each 
structure, and thus elucidated the mechanism of inducing the orbital polarization in the p symmetric 
states. We also discuss a simple sum rule connecting the XMCD spectra with the orbital moment 
in the p symmetric states. 

PACS numbers: 78.70.Dm, 71.20.Be, 75.50.Cc 



I. INTRODUCTION 

X-ray magnetic circular dichroism (XMCD) has at- 
tracted much interest as a useful tool to investigate mag- 
netic states!. The spectra directly reflect the magnetic 
order at the Z/2,3 edges of transition-metal compounds 
and at the edges of rare-earth compounds, be- 

cause by photoabsorption the core electron enters the 
3d or 4/ states which constitute the magnetic order. 
Sum rules have been found useful to evaluate the or- 
bital magnetization 2-3-4 . On the other hand, at earlier 
time, it was not clear where the XMCD spectra at the K 
edge of transition-metal compounds come from, because 
the Ap states the core electron enters by photoabsorption 
are not the states constituting the magnetic order^i 
One of the present authors and Hirai have analyzed the 
XMCD at the K edge of ferromagnetic metals, Fe, Co, 
and Ni, and have found that the spectra come from the 
Ap orbital polarization induced by the mixing to the 3d 
states at neighboring sites&&. Subsequent experiments 
and calculations reached the similar conclusion for sev- 
eral transition-metal compoundsiifliiii^iiiiiii^iSiil 

The above mechanism is a consequence of an extended 
character of Ap states, and is closely related to the mech- 
anism of the resonant x-ray scattering (RXS) at the K 
edges of transition metals in several compounds 18 - 19 . The 
RXS intensity on superlattice spots was considered to be 
a direct reflection of the orbital order in LaMnOj— , but 
subsequent studies based on band structure calculations 
have revealed that the RXS intensity arises mainly from 
the lattice distortion which modulates the Ap states in 
the intermediate state of the second order process of the 
RXS 21 i 22 i 23 . Recently, ab initio calculations of the RXS 
spectra have been carried out for YTi03 and YVO3, hav- 
ing clarified the important role of lattice distortions 4 ^. 

Despite such progress, the underlying mechanism of 
inducing Ap orbital polarizations is not fully understood. 
The XMCD experimental data recently obtained in the 
ferromagnetic phase of MnsGaGSSiS 7 . seems suitable to 
make clear the mechanism, because considerable signals 



have been found not only at the K edge of Mn but also at 
the K edge of Ga. In this paper, we analyze the XMCD 
spectra at the K edges in M^GaC through an ab initio 
calculation. We use the Korringa-Kohn-Rostoker (KKR) 
method within the muffin-tin (MT) approximation 3 * in 
the local density approximation (LDA) scheme, and take 
account of the spin-orbit interaction (SOI). We neglect 
the core-hole potential in the final state. We expect 
that the final-state interaction is unimportant except for 
the Fermi-edge singularity, because it is known from the 
study of the XMCD in Fe, Co, Ni, that the spectra are 
well reproduced without taking account of the final-state 
interactionSiL2iS. This may come from the fact that the 
Ap states have small amplitude inside the MT sphere and 
thus are little subject to the core hole potential. 

In Sec. II, we calculate the electronic structure in a fer- 
romagnetic phase of MnaGaC. In Sec. Ill, we present the 
formulas for the XMCD spectra and derive a sum rule. 
We also discuss the calculated spectra in comparison with 
the experiment. Section IV is devoted to concluding re- 
marks. 



II. ELECTRONIC STRUCTURE 

This material takes an "inverse" perovskite 
structures 8 ., as schematically shown in Fig. ^ It 
has an interesting magnetic property; the antifer- 
romagnetic phase in low temperatures turns into a 
ferromagnetic phase in high temperatures through a 
first-order transition at 168 We carry out the 

band calculation assuming a ferromagnetic phase with 
the magnetization fixed to the opposite of the [111] 
direction, although an antiferromagnetic phase may be 
stabler than the ferromagnetic phase. By fixing the 
magnetization direction, we neglect the small magnetic 
anisotropy. We define the density of states (DOS) inside 
the MT sphere at a Mn site and at a Ga site. On Mn 
sites, the DOS projected onto the d symmetric states is 
dominating near the Fermi level. The s and p symmetric 
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FIG. 1: Inverse perovskite-type crystal structure. 



spectra are defined by 

/ioH - [/*+(w) + /i_(u;)]/2, (3.2) 
fj, c (ui) = \fi + (w) - n-(uj)]. (3.3) 

The XMCD intensity arises from the p orbital polariza- 
tion in the unoccupied states, as is obvious from Eq. Q3.3JI . 
This quantity can be connected with the orbital moment 
in the p symmetric states as follows. Let the wave func- 
tion be expanded in the MT sphere as 



0n,k(r) 



E*n,k ,/;™' k (VI 



(3.4) 



DOS's are found extremely small. The spin and orbital 
angular momenta, denoted as S and L in units of h, are 
defined also inside the MT sphere. We get S — 0.73 and 
L = 0.042 in the d symmetric states, both of which are 
pointing to the [111] direction, while we get S — —0.0006, 
L = 0.00006 in the p symmetric states (the — sign means 
that the angular momentum is pointing to the opposite 
of the [111] direction). The total magnetic moment 
becomes —(25 + L)fiB = — 1.50/x.b per Mn. This value 
is consistent with the magnetization measurement^ and 
the previous band calculations without taking account 
of the SOI 2 ^. On the other hand, on Ga sites, the 3d 
states are located about 12 eV below the Fermi level. 
The s and p symmetric DOS's are found small above the 
Fermi level. Small angular momenta are induced due to 
the influence of the magnetic moments at Mn sites. We 
get S = —0.011, L = 0.00046 in the p symmetric states, 
and much smaller values, S = 0.0021, L = -0.00041, in 
the d symmetric states. 



III. XMCD SPECTRA 

We calculate the absorption coefficient, neglecting the 
core hole potential. Assuming that photons are propa- 
gating along the [111] direction, we have the expressions 
for the right-handed (+) and left-handed (— ) circular po- 
larizations, 



H±(w) <x J r2dr J dn<p n>k (r)*rY li±1 (n)R ls (r) 



x 8{fwj - e„,k)0(en,k - e F ), 



(3.1) 



where the spherical harmonic function Yi,±i is defined 
with the quantization axis along the [111] direction. The 
step function 9(x) ensures that the sum is taken over 
states above the Fermi level. The R\ s represents the Is 
wave function of Mn or Ga, and (f> n ,k represents the wave 
function with the band index n, wave- vector k and energy 
e n> k. In the actual calculation, we replace S(hu> — e n ,k) by 
a Lorentzian form (T / tt) / ((hu; — e„ ik ) 2 + T 2 ) with F = 1 
eV, in order to take account of the Is core hole life-time 
width. The total absorption coefficient and the XMCD 



where ip™' m CT ( r ) represents the wave function projected 
onto the state with the angular momentum {£,m}, and 
spin cr, and normalized inside the MT sphere. Then the 
orbital angular momentum in the p symmetric state is 
expressed as 



7i , k , a 



mi 



n,k |2 



l/£-l,J a )0(*F-*n,k), (3-5) 



where the sum is taken over the occupied states. Note 
that, if the sum is taken over all the states both oc- 
cupied and unoccupied, one may get Y^ n k <r(\fl'l a- 1 2 — 
\fi-i o-| 2 ) = 0. On the other hand, the absorption coef- 
ficient may be rewritten as 



/i±(w) oc ^ ~ en,k)0(en,k - ef)|/"±i 



Jr 2 dr f dn^ lt<T (r)*rY lt± i(Sl)R ls (r 



n,k |2 



.(3.6) 



The last factor represents the transition probability from 
the Is state to the p symmetric states, which varies lit- 
tle with respect to the energy of the p symmetric states 
above the Fermi level up to ~ 20 eV (within ~ 10% 
variation). This is because the p symmetric states are 
well approximated by the atomic 4p wave function on 
the region close to the center of atoms. Replacing the 
last factor by the value at the Fermi energy (denoted as 
|iSf| 2 ) in Ea. f3.fip . we obtain approximately the relation, 



fj, c (ui)dui oc — |S , f| 2 ^p 



(3.7) 



We need to introduce some upper cutoff a->2, because the 
transition probability is replaced by a constant |<Sf| 2 - It 
may be set by Tuo% = ep + 20 eV. This sum rule is an 
extension of the relation previously derived on the basis 
of the tight-binding model for Fe, Co, Ni metalsS^. 

Now we discuss the calculated spectra at the Mn K 
edge. Figure[2fa) and (b) show the total absorption coef- 
ficient ho(uj) and the XMCD spectra /i c (w) in comparison 
with the experiment? 2 -^. The experimental data are for a 
powdered sample under a small magnetic field H = 0.6 
T at 200 K. The sample is in the ferromagnetic phase 
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FIG. 2: (a)Total absorption coefficient jk>(w), and (b)XMCD 
spectra fj, c (u) at the Mn _K" edge, as a function of photon 
energy. The origin of energy corresponds to the excitation to 
the Fermi level. The normalized intensities are magnified to 
10 3 times. The experimental data are for a powdered sample 
in the ferromagnetic phase at T — 200 K 26 . (c)XMCD spectra 
calculated with turning off the SOI on all the states at Ga 
sites, on the p symmetric states at Mn sites, and on the d 
symmetric states at Mn sites, respectively. £a6 = means 
that the SOI on the b symmetric states at A sites is turned 
off. 



with the magnetization ~ 0.8/iB per Mn, which value is 
smaller than the calculated one 1.5/iB. The difference 
may be ascribed to a temperature effect, which usually 
reduces the magnetization. The calculated value of /x c (u;) 
is divided by the value of (Xq (w) at Hu = ep + 20 eV, while 
the experimental XMCD spectra are divided by the value 
of the total absorption coefficient at the energy about 40 
eV higher than the threshold. Both /zo(w) and ^ c (w) are 
in good agreement with the experiment. For the total ab- 
sorption coefficient, peaks A and B correspond well with 
shoulders A' and B' in the experimental curve. For the 
XMCD spectra, dips C, D, and peak E correspond well 
with the experimental ones C, D', and E'. Only one de- 
viation from the experiment is that the small peak at the 
Fermi level is not reproduced. The intensity integrated 
up to lu 2 is found negative. This is consistent with the 
sum rule, Ea. (|3.7|) . because the value of L p is positive at 
Mn sites. 



The XMCD spectra come from the orbital polarization 
in the p symmetric states, which may be induced by (i) 
the spin polarization in the p symmetric states through 
the SOI, and (ii) the orbital polarization at neighbor- 
ing sites through hybridization. For making clear the 
inducing mechanism, we calculate the XMCD spectra at 
Mn sites with turning off the SOI on several specified 
states. The top panel among three panels of Fig. J2Jc) 
shows the XMCD spectra with turning off the SOI on 
all the states at Ga sites. The spectra remain similar 
except for a slight suppression of peak E, indicating that 
the orbital polarization at Ga sites have little influence 
on the XMCD spectra at Mn sites. The middle panel 
in Fig. Efc) shows the XMCD spectra with turning off 
the SOI only on the p symmetric states at Mn sites. Dip 
C keeps the similar shape, while dip D almost vanishes. 
This indicates that the p orbital polarization correspond- 
ing to dip D is induced by the spin polarization in the 
p symmetric states through the SOI. The bottom panel 
shows the XMCD spectra with turning off the SOI only 
on the d symmetric states at Mn sites. The intensity of 
dip C is drastically reduced, but dip D and peak E re- 
mains similar, indicating that the 3d orbital polarization 
gives rise to dip C. Within the MT approximation, the 
3d orbital polarization cannot polarize the p orbital in 
the same Mn site, because the p-d Coulomb interaction 
is spherically averaged inside the MT sphere. Thus we 
conclude that the 3d orbital polarization at neighboring 
Mn sites induces the p orbital polarization corresponding 
to dip C through the p-d hybridization. 

Next we discuss the XMCD at the Ga K edge. Fig- 
ure [Sf a) and (b) shows the calculated /j,o(lj) and /x c (w) 
in comparison with the experimental. Both ho(oj) and 
fi c (w) are in good agreement with the experiment. For 
the total absorption coefficient, peaks F and G corre- 
spond well with the experimental ones. Also, for the 
XMCD spectra, dip H, peak I and dip J correspond well 
with the experimental ones. The intensity integrated up 
to L02 is found negative. This is again consistent with the 
sum rule, because the value of L p is positive at Ga sites. 

For making clear the mechanism, we also calculate the 
XMCD spectra at Ga sites with turning off the SOI on 
several specified states. The top panel in Fig.^c) shows 
the XMCD spectra with turning off the SOI on the p 
symmetric states at Ga sites. Peak I and dip J are al- 
most suppressed, indicating that the p orbital polariza- 
tion corresponding to those structures is induced by the 
spin polarization in the p symmetric states through the 
SOI. On the other hand, dip H is half suppressed, indi- 
cating that the above mechanism works in a small part 
on dip H. The bottom panel shows the XMCD spectra 
with turning off the SOI on the d symmetric states at Mn 
sites. Dip J remains similar, dip H is half suppressed, and 
peak I is even enhanced. This indicates that the nega- 
tive values are added to both dip H and peak I due to 
the coupling to the 3d orbital polarization at neighboring 
Mn sites. 
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FIG. 3: (a)Total absorption coefficient /uo(w),and (b)XMCD 
spectra /i c (u>) at the Ga K edge, as a function of photon en- 
ergy. The origin of energy corresponds to the excitation to 
the Fermi level. The normalized intensities are magnified to 
10 3 times. The experimental data are for a powdered sample 
in the ferromagnetic phase at T = 200 K 27 . (c)XMCD spec- 
tra calculated with turning off the SOI on the p symmetric 
states at Ga sites, and on the d symmetric states at Mn sites, 
respectively. £a& = means that the SOI on the b symmetric 
states at A sites is turned off. 



As shown above, the calculation reproduces well the 
experimental curves. Only the exception is the small 
positive peak W at the Fermi level on Mn sites in the 
XMCD spectra (see Fig. 2). We might take it a small er- 
ror, since the peak is small. However, in case of MnsZnC, 
the similar positive peak at the Fermi level becomes large 
on Mn sites. To gain confidence in our calculation, we 
carry out a similar calculation for M^ZnC, assuming a 
ferromagnetic phase. We get S — 1.088, L = 0.019 in 
the d symmetric states, and S = -0.0003, L = -0.00005 
in the p symmetric states on Mn site. Figure 0] shows 
the calculated absorption and XMCD spectra at Mn K 
edge, in comparison with the experiment at T = 300 K 
(ferromagnetic phase). The positive peak at the Fermi 
level is clearly reproduced. 



IV. CONCLUDING REMARKS 

We have reported an ab initio calculation of the XMCD 
spectra at the K edges of Mn and Ga in the ferromag- 
netic phase of M^GaC, by taking account of the SOI 
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FIG. 4: (a)Total absorption coefficient /xo(w),and (b)XMCD 
spectra /U c (w) at the Mn K edge in MnaZnC, as a function 
of photon energy. The origin of energy corresponds to the 
excitation to the Fermi level. The normalized intensities are 
magnified to 10 3 times. The experimental data are for a pow- 
dered sample in the ferromagnetic phase at T = 300 K^. 



in the LDA scheme. The calculated spectra show excel- 
lent agreement with the recent experiment. The spec- 
tra have explicitly been shown to arise from the p or- 
bital polarization, and the mechanism of its induction 
has been fully analyzed; the associated process is identi- 
fied for each structure. The present result may serve a 
guide to analyze the K edge XMCD spectra in other ma- 
terials. Such usefulness has been demonstrated for the 
ferromagnetic Ni metal. We have also derived a simple 
sum rule, which connects the XMCD spectra with the 
orbital moment, and have confirmed it works. Judging 
from the agreement of the calculated spectra with the 
experiment, we think that the effect of the ls-core-hole 
potential in the final state is negligibly small. 

The present finding is rather general and can be ap- 
plied to other systems. To demonstrate this, we compare 
the XMCD spectra in M^GaC with those in the ferro- 
magnetic Ni metal. Figure |S] shows the latter quantity 
calculated with the same method as above. We divide 
/i c (w) by iM)(u>2) at T1UJ2 = ef + 20 eV. A large dip K 
appears above the Fermi level; it is half suppressed if the 
SOI is turned off either in the p symmetric states or in the 
3c? states. This behavior is similar to dip C in MnaGaC, 
although the contribution of the SOI in the p symmetric 
states is a little larger in Ni. On the other hand, a tiny 
dip L at high energy disappears without the SOI in the 
p symmetric states, similar to dip D in M^GaC. 

Finally we point out that this type of ab initio calcula- 
tion can be applied to analyzing the magnetic RXS spec- 
tra at the K edge in transition-metal compounds^iiS, 
and may clarify the p orbital polarization in antiferro- 
magnetic phases. According to our recent calculation of 
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the magnetic RXS spectra in NiO^, the main peak in- 
tensity almost disappears when the SOI on the 4p band 
is turned off. This finding is consistent with the above 
analysis, because the 4p states are located at ~ 10 eV 
higher than the Fermi level. 
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